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Abstract 

The interface problem describing the scattering of time-harmonic electromagnetic 
waves by a dielectric body is often formulated as a pair of coupled boundary integral 
equations for the electric and magnetic current densities on the interface T. In this 
paper, following an idea developed by Kleinman and Martin [H] for acoustic scatter- 
ing problems, we consider methods for solving the dielectric scattering problem using 
a single integral equation over T for a single unknown density. One knows that such 
boundary integral formulations of the Maxwell equations are not uniquely solvable 
when the exterior wave number is an eigenvalue of an associated interior Maxwell 
boundary value problem. We obtain four different families of integral equations for 
which we can show that by choosing some parameters in an appropriate way, they 
become uniquely solvable for all real frequencies. We analyze the well-posedness 
of the integral equations in the space of finite energy on smooth and non-smooth 
boundaries. 

Keywords : scattering problems, Maxwell equations, boundary integral equations, 
Helmholtz decomposition. 

1 Introduction 

We consider the scattering of time-harmonic electromagnetic waves in M 3 by a bounded 
Lipschitz obstacle. We assume that the dielectric permittivity and the magnetic perme- 
ability take constant, in general different, values in the interior and in the exterior of 
the domain. This problem is described by the system of Maxwell's equations, valid in 
the sense of distributions in R 3 , which implies two transmission conditions expressing the 
continuity of the tangential components of the fields across the interface. The transmis- 
sion problem is completed by the Silver-Miiller radiation condition at infinity (see [26] 
and [28]). 

It is well known that this problem can be reduced in several different ways to systems 
of two boundary integral equations for two unknown tangential vector fields on the inter- 
face. Such formulations are analyzed in Harrington's book [16] and in Martin and Ola's 
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comprehensive paper |22j. Some pairs of boundary integral equations, such as Midler's 
|27| . are uniquely solvable for all real values of the exterior wave number and others, such 
as the so-called electric- field formulations [22j are not, although the underlying Maxwell 
interface problem is always uniquely solvable under standard assumptions on the material 
coefficients. 

More recent research works in the scientific and engineering community show that 
there are computational advantages to solve dielectric scattering problems via a single 
integral equation for a single unknown, rather than a system of two equations of two 
unknowns. For two-dimensional dielectric scattering problems, one can find various for- 
mulations and numerical results in [29, 30, 32j. In |23j Marx develops single source integral 
formulations for three-dimensional homogeneous dielectric objects using an ansatz on the 
exterior electric field and in |34j Yeung presents electric-field (EFIE) and magnetic-field 
(MFIE) integral equations for a single unknown based on an ansatz on the interior scat- 
tered field. Computational results in [34j show higher convergence speed for the MFIE 
and the EFIE than for the pairs of boundary integral equations. However, both of these 
single integral formulations suffer from spurious non-unique solvability due to interior 
resonances. 

In this paper, we study methods for solving the transmission problem using a single 
boundary integral equation for a single unknown tangential vector field on the interface 
by eliminating irregular frequencies. We follow ideas of |19j where Kleinman and Martin 
considered the analogous question for the acoustic interface scattering problem. The 
method consists of representing the solution in one domain by some combination of a 
single layer potential and a double layer potential, and inserting this representation into 
the transmission conditions and the Calderon relations of the other domain. Several 
different integral equations of the first kind or of the second kind, containing two arbitrary 
parameters, can be obtained in this way, and in the scalar case, the parameters can be 
chosen in such a way that no spurious real frequencies are introduced. Following the same 
procedure in the electromagnetic case, one encounters two main difficulties: 

The first problem is that some boundary integral operators that are compact in the 
scalar case are no longer compact, and therefore arguments based on the theory of Fred- 
holm integral equations of the second kind have to be refined in order to show well- 
posedness of the corresponding integral equations. 

The second problem comes from a lack of ellipticity. The spurious frequencies are 
associated with the spectrum of a certain interior boundary value problem of the third 
kind, and whereas in the scalar case this is an elliptic boundary value problem whose 
spectrum can be moved off the real line by the right choice of parameters, in the Maxwell 
case this boundary value problem is not elliptic, in general. Thus an additional idea is 
needed to avoid real irregular frequencies. Mautz presents in [2l| an alternative that leads 
to an associated interior problem with an impedance boundary condition, which ensures 
the uniqueness of the solution. However Mautz's equation is not adapted to our point of 
view of variational methods and energy spaces. Since the Kleinman-Martin method has 
similarities to the combined field integral equation method, we use a regularizer introduced 
by Steinbach and Windisch in [31] in the context of combined field integral equations for 
the time-harmonic Maxwell equations. This regularizer is a positive definite boundary 
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integral operator with a similar structure as the operator of the electrical field integral 
equation, but it is not a compact operator like those used in |10j and [6] for regularizing 
the exterior electromagnetic scattering problem. Its introduction changes the boundary 
condition in the associated interior boundary value problem from a non-elliptic local 
impedance-like condition to a non-local, but elliptic, boundary condition. 

This work contains results from the thesis [21J where this integral formulation of the 
transmission problem is used to study the shape derivatives of the solution of the dielectric 
scattering problem, in the context of a problem of optimizing the shape of a dielectric 
lens in order to obtain a prescribed radiation pattern. 

In Section [3] we recall some results about traces and potentials for Maxwell's equations 
in Sobolev spaces. We use the notation of [7] and |6j and quote some important properties 
of the boundary integral operators that constitute the Calderon projector for Maxwell's 
equations. 

Sections H] and [5] contain the details of the method for solving the transmission prob- 
lem using single-source boundary integral equations. In Section [H we start from a layer 
representation for the exterior field whereas in Section O we use a layer representation 
for the interior field. In either case, we derive two boundary integral equations of the 
second kind and we show uniqueness of their solutions under suitable conditions on an 
associated interior boundary value problem. Moreover, we show that the integral opera- 
tors in each integral equation are Fredholm of index zero. We also construct the solution 
of the transmission problem using the solution of any of the four integral equations. We 
finally show how to choose the free parameters so that the associated interior boundary 
value problem is uniquely solvable, and as a consequence, we can construct an integral 
representation of the solution which yields uniquely solvable boundary integral equations 
for all real frequencies. 

For smooth domains, we base the analysis of the integral operators on the technique 
of Helmholtz decomposition, which represents a tangential vector field by two scalar 
field and each integral operator acting on tangential fields by a two-by-two matrix of 
scalar operators. Since these operators then act between standard Sobolev spaces instead 
of the complicated mixed-order energy space, it is easy to check for compactness or 
ellipticity. Using this technique, we find rather general sufficient conditions on the physical 
parameters to ensure unique solvability of the integral equations. If the boundary is only 
Lipschitz, we show that under more restrictive conditions one still has strong ellipticity 
of the integral operators. This conditions include the physically relevant case of positive 
permeabilities, permittivities, and frequencies. 

2 The dielectric scattering problem 

Let O denote a bounded domain in M 3 and let f2 c denote the exterior domain M 3 \fi. In 
this paper, we will assume that the boundary r of fHs a Lipschitz continuous and simply 
connected closed surface. Let n denote the outer unit normal vector on the boundary T. 

In f2 (resp. £l c ) the electric permittivity e« (resp. e e ) and the magnetic permeability 
/ij (resp. (i e ) are positive constants. The frequency u is the same in f2 and in f2 c . The 
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interior wave number m and the exterior wave number n e are complex constants of non 
negative imaginary part. 

Notation: For a domain G C t 3 we denote by H S (G) the usual L 2 -based Sobolev 
space of order s G R, and by Hf QC (G) the space of functions whose restrictions to any 
bounded subdomain B of G belong to H S (B), with the convention H° = L 2 . Spaces of 
vector functions will be denoted by boldface letters, thus 

H S (G) = (H s (G)f. 

If D is a differential operator, we write: 

H(D, fl) = {u G L 2 (n) : Du G L 2 (f])} 
H loc (D,JF) = {u G Lf oc (?F) : Dn G L 2 oc (W)} 

The space H(D,il) is endowed with the natural graph norm. This defines in particular 
the Hilbert spaces H(curl, $7) and H (curl curl, Q). 

The time-harmonic dielectric scattering problem is formulated as follows. 



The dielectric scattering problem : 

Given an incident_field E inc G H loc (curl, R 3 ) that satisfies curl curl E inc - K |E inc 
neighborhood of CI, we seek two fields E l G H(curl, ft,) and E s G Hi oc (curl, Q c ) satisfying 
the time-harmonic Maxwell equations 



K 2 E mc = 0ina 



curlcurir-K 2 ^ =0 







curl curl E s 

the two transmission conditions, 

n x V = n x (E s + E inc ) 
/tx~ 1 (n x curlE*) = u f T 1 n x curl(E s + E 

and the Silver-Muller radiation condition: 



in Q, 
in fi c , 



on r 

on r 
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curlE s (x) x — -in e E s (x) 



0. 



(2.1) 
(2.2) 



(2.3) 
(2.4) 



(2.5) 



It is well known that this problem has at most one solution under some mild restric- 
tions on the dielectric constants. We give sufficient conditions in the next theorem, and 
for completeness we give its simple proof. 

Theorem 2.1 Assume that the constants fj,i, Ki, fi e and n e satisfy: 
(i) n e is real and positive or Im(re e ) > 0, 



fii) Im ( K e — ] < and Im ( K e — K 2 ) > 0. 
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Then the dielectric scattering problem has at most one solution. 



PROOF. We use similar arguments as in the acoustic case |19j . Assume that E mc = 0. 
Let (E l , E s ) be a solution of the homogeneous scattering problem. Let B R be a ball of 
radius R large enough such that Vt C B R and let n R the unit outer normal vector to Br. 
Integration by parts using the Maxwell equations (|2.ip and (|2.2p and the transmission 
conditions fj2.3|) and f|2.4[) gives : 

I (curl E s x n R ) ■ W = [ {| curl E s | 2 - k 2 |E s | 2 } + ^ [ {| curl E*| 2 - k 2 )^ 2 } 

JdB R J Hi Jfl 

B R \n 

We multiply this by T^T and take the imaginary part: 

Imp^/ (curl E s x n R ) ■ W J = Im(T^) f f {| curl E s | 2 + \n e E s \ 2 } ) 
\ JdB R J \JB R \n J 

+ l m (-A) [ |curlEf -Imf-A* 2 ) [ |Ef. 
V IH J Jn \ ih J J n 

Under the hypotheses (i) and (ii), all terms on the right hand side are non-positive. 
Thanks to the Silver-Muller condition, we have 

lim / | curl E s x n R - in e E s \ 2 = 0. 

R^+oo Jq Br 

Developing this expression, we get 



lim 



ldB R 

As we have seen, we have 



/ |curlE s x n fi | 2 + |«: e E s | 2 -2 Re (curl E s x n R ■ w e E s ) = 0. 



/ Re (curl E s x ■ in e E) = Im / (n e curl E s x riR • E) < 0. 



>dB R JdB R 
It follows that 

lim / |E S | 2 = 0. 



R-^+oo Jq Br 

Thus, by Rellich's lemma |10| . E s = in Sl c . Using the transmission conditions, we obtain 
lD E { = j Nk . E i = 0. It follows that E* = in 9,. U 



3 Traces and electromagnetic potentials 

We use some well known results about traces of vector fields and integral representations 
of time-harmonic electromagnetic fields on a bounded domain Q. Details can be found in 
HI El [3 [TTJ [28] . Recall that the boundary T is only assumed to be Lipschitz continuous, 
unless stated otherwise. 



5 



Definition 3.1 For a vector function u € (^°°(f2)) 3 and a scalar function v G < ^' 00 (f2) 
we define the traces : 

•yv = v\ r 
715U = (n x u)| r (Dirichlet) 
Jn k u = k~ (n x curlu)| r (Neumann). 

We use standard Sobolev spaces if*(T), £ € [—1,1], endowed with standard norms 
II ' ll.ff*(r) an d with the convention H°(T) = L 2 (F). Spaces of vector densities are denoted 
by boldface letters, thus H*(r) = (H^T)) . We define spaces of tangential vector fields 
as H x (r) = n x H*(T). Note that on non-smooth boundaries, the latter space different, in 
general, from the other space of tangential vector fields Hjj(r) = n x H* X (T). On smooth 
boundaries, the two spaces coincide. For s = we set H x (r) = L x (r). 

The trace maps 

7 : H s+ ^(n) -> H S (T), 

are continuous for all s > 0, if the domain is smooth. On a polyhedron, the trace maps 
are continuous for s £ (0, 2), whereas for a general bounded Lipschitz domain in general, 
the validity is only given for s € (0, 1). For s = 1, the trace operator 7 fails, in general, 
to map Hz(Q) to -ff 1 (r), although H 1 ^) is well defined on the boundary T, see [18J. 

The dual spaces of i?*(r) and H* x (r) with respect to the L? (or L 2 ) scalar product is 
denoted by H~ t (F) and H~'(r), respectively. 

We use the surface differential operators: The tangential gradient denoted by Vr, the 
surface divergence denoted by divr, the tangential vector curl denoted by curlr and the 
surface scalar curl denoted by curlr- For their definitions we refer to [5], |11| and [28J. 

Definition 3.2 We define the Hilbert space 

H7(div r ,r) = jj G H7(r),div r j e ff-s(r) 

endowed with the norm 

1 1 • 1 1 1 = 1 1 ■ 1 1 1 +11 divr • 1 1 1 

H x ?(div r ,r) H x ?(r) n H-?(F) 

The skew-symmetric bilinear form 

B: (div r , T) x (div r , T) -> C 

( j, m) — y B( j, m) = J j -(m x n) da 

_i 

defines a non-degenerate duality product on H x 2 (divr,r). 



6 



Lemma 3.3 The operators jrj and are linear and continuous from (^°°(J7)) 3 to 

Lx(r) and they can be extended to continuous linear operators from H(curl,il) and 

_i 

H(curl, O) fl H (curl curl, fl), respectively, to H x 2 (divr,r). Moreover, for all u, v G 
H(curl, 0), we /iat>e: 



/ 



[(curl u • v) — (u • curl v)] dx = B{^dv, 7du). (3.1) 



For u G Hi oc (curl, f2 c ) and v € Hi oc (curl curl, S7 C )) we define 7^,11 and ^ c N \i in the 
same way and the same mapping properties hold true. 

Let k be a complex number such that Im(«;) > and let 

gin\x— y\ 

G{K,\x-y\) = — 1 

4tt\x — y\ 

be the fundamental solution of the Helmholtz equation 

An + k 2 u = 0. 
The single layer potential ip K is given by : 

(il> K u)(x) = [ G{k, \x - y\)u(y)da(y) x G R 3 \T, 



and its trace by 

V R u{x) = I G(k, \x - y\)u(y)da(y) x G T. 



For a proof of the following well-known result, see \17\ 128] . 
Lemma 3.4 The operators 

V K : H-^T) ^ h'HF) 

are continuous. 

_ i 

We define the electric potential ^ e k generated by j G H x 2 (divr,r) by 

This can be written as ^e k j := K_1 curl curl ip K ] because of the Helmholtz equation and 
the identity curl curl = — A + Vdiv (cf. [3]). 

_ i 

We define the magnetic potential $m« generated by m G H x 2 (divr,r) by 

^j\/ K m := curiam. 

These potentials satisfy 

K curl^^ = ^m k and k~ curl^>M K = ^e k - 
We denote the identity operator by I. 
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Lemma 3.5 The potential operators ^e k an d^M K are continuous from H x 2 (divr, T) to 

_ i 

Hi oc (curl,R 3 ). For] G H x 2 (div r ,r) we have 

(curl curl -k 2 1)^e k \ = and (curl curl -k 2 I)^m k iti = in R 3 \r 
and ^ e k j an d ^M K m satisfy the Silver-Muller condition. 

It follows that the traces 70, "fN K7 7^ and 7^ can be applied to \& e k an d ^m^j resulting 

_ 1 

in continuous mappings from H x 2 (divr,r) to itself satisfying 

1N K ^E K = 1D^M K and 7AT K *M K = 1D^E K ■ 

Defining 

[7o] = 1D-1 C D , {ID} = -^(lD+l C D ), 

we have the following jump relations (see |7j): 

[7d]^s k = 0, [7jvJ = -I, 

[7d]*m k = -I, [7JV„] *Af„ = 0. 

Now assume that E € L^ oc (R 3 ) belongs to H(curl, f2) in the interior domain and to 
Hi oc (curl, f2 c ) in the exterior domain and satisfies the equation 

(curl curl -k 2 I) E = (3.2) 

in R 3 \ T and the Silver-Muller condition. Then if we set j = [7jvJE, m = [td]E, we have 
on R 3 \ T the Stratton-Chu integral representation 

E = --*eJ-*m k ™. (3.3) 
Special cases of (|3.3p are: If (E*, E s ) solves the dielectric scattering problem, then 

-^ Nmm E'-9 MHmih E' = [- E ' H% (3.4) 
*^ e 7k e ( ES + ^ + *M Ke lh (E s + P-) = I E £mc I I % (3.5) 

-*B Hi TN ni #-*M„ i lD# = [ % (3.6) 
We can now define the main boundary integral operators: 

C K = {7d}^e k = {7ivJ*M K , 
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M K = {7Z?}*M K = {lN K }^E K - 



_ l 

These are bounded operators in H x 2 (divr,T). 

As tools, we will need variants of these operators: 



Definition 3.6 Define the operators Mq, C Kj o and Cq for} € H x 2 (divr,r) by : 

M } = -{7 D }curl^ , 

C K ,o j = -« n x V j +KT 1 curl r V div r j, 
Cq j = n x Vq j + curlr Vb divr j . 

Note that Cq differs from Ci ; o by the relative sign of the two terms. 

We collect now some properties of these boundary integral operators that are known 
for Lipschitz domains. 

First we note the following useful relations: 

curlr Vr = and divr curlr = (3.7) 

divr(n x j) = — curlr j and curlr (n x j) = divr j (3.8) 
The following lemma is proved in [U [7] . 

_ i 

Lemma 3.7 (i) The operators C K — C Kj o and M K — Mq are compact in H x 2 (divr, T). 
(ii) Both C K and M K are antisymmetric with respect to the bilinear form B. 

The Calderon projectors for the time-harmonic Maxwell system ([3.2p are P = ^1 + A K 
and P c = f I — A.k where 

A= ( M R C R 
\C K M K 

We have P o P c = and therefore 

C 2 K = \l- M 2 K and C K M K = -M K C K . (3.9) 

It is a classical result that when the boundary V is smooth, the operator M K is compact 

_ i 

from H x 2 (divr,r) to itself (see [28J). From identity (|3.9p one can then immediately 
deduce that the "electric" operator C K is Fredholm of index zero. The latter result is 
also true for a Lipschitz boundary V (see (Tj for more details). Later on, we need the 
corresponding result for the "magnetic" operator ^1 + M K . This has been proved for 
Lipschitz domains in [25, Thm. 4.8] and in |31| Thm. 3.2]: 

_i _ i 

Lemma 3.8 The operator |l + M K : H x 2 (divr,r) — > H x 2 (divr,r) is Fredholm of index 
zero. 
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In fact, for IrriK > this operator is an isomorphism. Later on, we use the result for 
K = 0, where we don't know whether it is an isomorphism. But we only need that is an 
isomorphism up to a compact perturbation, that is, Fredholm of index zero, so we will 
not pursue this further. 

The following theorem was proved in |31| Thm. 2.6]. 



Lemma 3.9 The operator Cq is self-adjoint and elliptic for the bilinear form B and 
_i 

invertible on H x 2 (divr,r). Ellipticity means here that there exists a positive constant a 
_ i 

such that for all] G H x 2 (divp, T) 

B(C *j,j) >a||j|| 2 i 

H x 2 (div r ,r) 

_ l 

Indeed, for j € H x 2 (divr, T) we have 



BQ, Cfj) = J {j -Vq} + div r j Vb div r 1} (3.10) 

and the result follows from the H~2 (r)-ellipticity of the scalar single layer potential 
operator Vo- 

4 Integral equations 1 

In this section, we present the first method for solving the dielectric problem, following 
the procedure of R. E. Kleinman and P. A. Martin |19| : We use a layer ansatz on the 
exterior field to construct two alternative boundary integral equations. 
In the scalar case, one represents the exterior field as a linear combination of a single 
layer potential and a double layer potential, both generated by the same density. It turns 
out that this simple idea does not suffice in the electromagnetic case if one wants to avoid 
irregular frequencies. Our approach is related to the idea of "modified combined field 
integral equations": We compose one of the electromagnetic potential operators with an 
elliptic and invertible boundary integral operator, namely Cq. More precisely, we assume 
that E s admits the following integral representation : 

E s (x) = -a(^ EK J)(x) - b(? MKe C$ j)(x) for x G W . (4.1) 



_ i 

Here j G H x 2 (divr,r) is the unknown density and a and b are arbitrary complex con- 
stants. 

We set p = ^ e % . The transmission conditions can be rewritten : 

7^E l = 7f,E* + 7£) E inc and >y N E* = p~ x ( 7 C N E s + 7A r Ke E* 



Using this in the integral representation formula (|3.6p in f2, we get: 

E 4 = --*E Ki (lk e E S + 7N Ke E m j - (7dE s + 7 dE- c ) in O. (4.2) 
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We take traces in (|4.ip and obtain the Calderon relations 

7 f,E s = {aC Ke -b(\l-M Ke )C%}\ =L e jonr, (4.3) 
7 ^E S = {- a (\l-M Ke )+bC Ke C* }] ^iV e jonr. (4.4) 

On the other hand, taking traces in (I4.2|) gives: 

pf-ll + M^ {%E S + lD E inc ) + C Ki ^% k E s + lNn E mc ) = Q onT, (4.5) 

f-il + M Ki ) ( 7 ^ E s + 7JVKe E- c ) + ( 7 ^E S + 7D E mc ) = on T. (4.6) 

We can now substitute (|4.3p and (|4.4jl into (|4.5p and get our first integral equation: 

Sj = /£>r~I + M« i U e j+C f(i JV e j = / (4.7) 

where / = -p (~1 + M K )j lD E inc - C k . 1Nkb E inc . (4.8) 

If we substitute (|4.3|) and (|4.4p into (|4.6p . we get our second integral equation: 

T] = pC Ki L e ]+(- l -l + M K )jN e \ = g (4.9) 

where g = -pC KilD E inc - (-±1 + M K )j 7J y Ke E mc . (4.10) 

Thus we obtain two boundary integral equations for the unknown j. Having solved 
either one, we construct E s using (|4.ip and E* using (I4.2p . f|4.3[) . (|4.4p : 

B = ~\ ^E Ki {lN K E™ + N e \}) _ (^{ 7D E mc + L e j}) . (4.11) 

_ i 

Theorem 4.1 If j G H x 2 (div r ,r) solves (|4~7) or (f4~9|) . i/ien E s and E* given by (|4~T]l 
cme? ()4.1ip solve the transmission problem. 

PROOF. We know that E l and E s satisfy the Maxwell equations and the Silver-Miiller 
condition. It remains to verify that E s and E* satisfy the transmission conditions (|2.3p and 
(|2,4p . Using the integral representation (|4.ip and (|4.2p of E s and E*, a simple computation 
gives: 

p( 1 %E s + lD E inc - lD E i ) = S\-f (4.12) 

and 

lN Ke E S + lN Ke E mc - pj NKi e = T j -g (4.13) 

We deduce that 

- if j solves (|4.7|) . then relation ([4.12p proves that the condition (|2.3p is satisfied, 

- if j solves (|4.9p . then relation (|4.13p proves that the condition (|2.4p is satisfied. 

Now we show that (|4.7p and (|4.9p are in fact equivalent. Define : 
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u(x) = { 7 iV Ke e nc + N e j}(x) - p* Mk . {7DE inc + L e j}(x) for x G ft c . 

This field u is in Hi oc (curl, f2 c ) and satisfies the Maxwell system 

curl curl u — nfu = (4-14) 
in Q c . On the boundary V we have: 

7o" = Sj-/ and 1n k . u = Tj —g. 

t 

Since u solves (I4.14p in Q c and satisfies the Silver-Miiller condition, it follows: 

j satisfies (|4.7p 7^11 = u = in 7 ^ u = =H satisfies (|4.9jl . 

j satisfies (f4~9|) 7^ u = => u = in tt c j^u = => j satisfies (|4.7|) . 

As a consequence, if j solves one of the two integral equations, it solves both, and then 
both transmission conditions (12.31) and (12.41) are satisfied. ■ 



The next theorem is concerned with the uniqueness of the solutions of the boundary 
integral equations (j4.7|) and (I4.9p . i.e., with the existence of nontrivial solutions of the 
following homogeneous forms of (|4.7p and fj4.9j> : 

p + M Ki ) L e j +C Kl iV e j = 0, (4.15) 
P C Ki L e j + f-ll + iV e j = 0. (4.16) 



As in the scalar case |19j . we associate with the dielectric scattering problem a new 
interior boundary value problem, the eigenvalues of which determine uniqueness for the 
integral equations. 

Associated interior problem: For a, b £ C, consider the boundary value problem 

curl curl u — KgU = in f2, cljdu — 6Co 7 Ar e u = on T. (4.17) 
Lemma 4.2 . Let a, b € C \ {0} and let K e £ C. Assume that 

(l) Im(^) ^0 i/Ke^K, 



(Mp Im(/tg) • Im (^ K e-^J > if K e £ 
then Kg is not an eigenvalue of the the associated interior problem (|4.17p . 

PROOF. Let Kg be an eigenvalue of the interior problem and let u 7^ be an eigenfunc- 
tion. Using Green's theorem we have: 

/ |curlu| 2 -Kg/ |u| 2 = K e ^( 7 DU,7Af Ke U) = K e S(7 jD U, 7 Ar Ke U) 

J(] Jn 

= (k 6 ^) S( 7d u, (Cq*)- 1 ^) if b + 

= 7^~B(Coj N u, 7jv Ke u) if a / 
a 
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Since Cq is elliptic for the bilinear form B, taking the imaginary part, we obtain 
-Im( K 2)^| u |2 =-Im( Ke ^)^((C *)- 1 7DU,7^u) 

= -|K e | 2 Im ( ^- ] S(7 JVfte u,C *7^-u). 



K e a t 

Under the hypotheses of the lemma the left-hand side and the right-hand side have op- 
posite sign, and it follows 

^((^o) _1 (7du),7^u) = and B(j Nk u, C *7aF^u) = 0. 

As Cq is elliptic for the bilinear form B, the traces 7dU and 7at K( ,u then vanish. Thanks 
to the Stratton-Chu representation formula (|3,3p in ft, we deduce that u = 0, which 
contradicts the initial assumption. ■ 

Remark 4.3 Note that this associated interior problem is not an impedance problem (or 

Robin problem) as in the scalar case flty . If we replace in (|4,17p the operator Cq by the 

identity, we obtain a "pseudo-impedance" type problem. This is a non-elliptic problem, 

about whose spectrum we have no information. That the problem is non-elliptic can be 

seen as follows: If it were elliptic, its principal part would be elliptic, too. This would be 

the vector Laplace operator with the "Neumann" condition 77v«. e U = 0. Any gradient of a 

harmonic function in H 1 (ft) will satisfy the homogeneous problem, which therefore has an 

infinite- dimensional nullspace, contradicting ellipticity. Note that the issue here is not the 

apparent non-elliptic nature of the interior Maxwell curl curl operator, which can easily 

be remedied by the usual regularization that adds — V div, but the manifestly non-elliptic 

nature of the Maxwell "Neumann" boundary operator. For a "true" impedance problem, 

the operator Cq would have to be replaced not by the identity, but by the rotation operator 

j i— y n X j which is used in Mautz's formulation WM- This operator leads out of the space 
_ i 

H x 2 (divr,r), however, which rules it out for our purposes. 



For our integral equations, the problem (|4,17p plays the same role as the Robin prob- 
lem for the scalar case in |19|. 



Theorem 4.4 Assume that the hypotheses of Theorem \2.1\ are satisfied. Then for (a, b) ^ 
(0,0), the homogeneous integral equations (|4.15p and (|4.16p admit nontrivial solutions if 
and only if is an eigenvalue of the associated interior problem. 

PROOF. Assume that j / solves (|4.15p or (|4.16p . 
We construct 112 and Ui as follows: 

u 2 (x) = -a^ EK J (x) - b^M Ke C$} (x) for x £ Q c 

m (a?) = (iV e j o )0c) - ^M Ki {L e ) Q )(x) for x € ft 

By Theorem 14. 1| Ui and 112 together solve the transmission problem with E mc = 0. 
Since this problem admits at most one solution, we have 112 = in Q c and Ui = in ft. 



13 



Now we set u(x) 
We have on T : 



a®E Ke io(x) - b^M Ke C^} (x) for x € Q 



7d u 2 -Idu = bC^ j , 
l C N Ke »2 -lN Ke u = aj . 



(4.18) 
(4.19) 



Since j c D u 2 = 7Ar e u 2 



on r, we find 



ajr,u — 6Co7Ar Ke u = on 17 



Thus u is an eigenfunction associated with the eigenvalue k 2 of the interior problem or 
u = 0. But this latter possibility can be eliminated since it implies that = Jn Kb u = 0, 
whence j = by (|4,18p and (I4.19p . which is contrary to the assumption. 
Conversely, assume that Kg is an eigenvalue of the associated interior problem. Let Vo ^ 
be a corresponding eigenfunction. The Calderon relations in f2 imply that : 



Using the equality a^ovo — b Co7./v Ke Vo = 0, we obtain 

L e (Co) _1 7z>wo = 0, N e {C* )- 1 lD v Q = 0, L e -y NKe v o = 0, N elNKe v = 0. 

If b 7^ 0, then 7_Dt>o 7^ 0, and j = (Cq) -1 7dUo is a nontrivial solution of (|4.15|) and (|4.16p . 
If b = 0, then ~fN Ke Vo 7^ 0, and j = ~fN Ke Vo is a nontrivial solution of (|4.15p and (|4.16p . ■ 

Whereas until now we only assumed that the boundary T is Lipschitz, we now present 
two theorems on the operators S and T that are, in this generality, only valid for smooth 
boundaries. By smooth we mean, for simplicity, < ^°° regularity, although a careful check 
of the proof would show that some finite regularity, such as c tf 2 , would be sufficient. 

Theorem 4.5 Assume that 

(i) the boundary T is smooth and simply connected, 
(ii) the constants a, b, p, e , /Uj, n e and m satisfy: 



-C Ke -/N K „ v + (5I - M Ke ) 7£)U = 
(|l - M Ke ) 7iv Ke v - C K ^ D v Q = 0. 




) 



^0. 



Then S is a Fredholm operator of index zero on H x 2 (divr,r). 



PROOF. We can rewrite S as follows: 



S = \bpC* Q - \bp{M K% + M Ke )C * + bpM Ki M Ke Cq - \ap{C Ke - C Ke , ) 
-\a{C Ki - C Kifi ) - \a{pC Kafi + C Kifi ) + apM Ki C Ke + aC Ki M Ke 
+b{C Ki — C Kii o)C Ke Cg + bC Ki: o(C Ke — C Kei o)Co + bC Ki fiC Kei oCQ. 
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Thus S is a compact perturbation of the operator 

Si = b {\pl + C Kifi C Kefi ) C * - -a (pC Kefi + C Kii0 ) ■ 

We have to show that the operator Si is Fredholm of index zero. For this we use the 

_ i 

Helmholtz decomposition of H x 2 (divr,r) : 

H~* (div r , T) = V t hI (T) curlr (r) . (4.20) 

For a detailed proof of (I4.20P see |15j . Note that we are assuming that the boundary T is 
smooth and simply connected. For a proof of the following result, we refer to [3] UH [28], 

Lemma 4.6 Let T be smooth and simply connected and t € R. The Laplace- Beltrami 
operator 

Ap = divr Vr = — curlr curlr (4-21) 

is linear and continuous from H t+2 (T) to i?*(r). 

It is an isomorphism from ilf* +2 (r)/R to the space Hl(F) defined by 

ueHt(T) ueH\T)and J u = 0. 



The terms in the decomposition j = Vrp + curlr Q f° r j 6 H x 2 (divr, T) are obtained by 
solving the Laplace-Beltrami equation: 



p = Ap 1 divr j , q = — Ap 1 curlr j . 

The mapping 



H x 2 (div r , T) -> H% (V) /M x (T) 

p \ (4-22) 



j = VyP + curlr q h ^ 

is an isomorphism. Using this isomorphism, we can rewrite the operator Si as an operator 
iSi defined from H a(r)/lR x H2(T)/M. into itself. Then to show that Si it is Fredholm of 
index zero it suffices to show that S± has this property. Let us begin by rewriting Cq and 

3 1 

C K; o as operators Cq and C Ki o defined on #2(r)/IR x f/~2(r)/IR. We have to determine 

P € H*(T)/R and Q G H?(T)/R such that C$(V r p + curl r <?) = Vr^o + curl r Qo, 
and this defines Cq by: 

(p \ _ f Po 

We have 



Ap 1 div r Cq(VtP + curl r < 



and 

Qo = -ApT 1 curlr C *(V r p + curl r q). 
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Using the integral representation of Cq and the equalities (|3.7p and fj3.8j) we obtain: 



r 



° V ^21,1 + ^21,2 ^22 

where 

Cn = —Ap 1 curlr VbVr, C12 = —Ap 1 curlr Vb curlr, 
C 2 i,i = -Ap 1 div r V Vr, C22 = -Ap 1 div r V" curl r , 
C2i,2 = UoAr. 

Some of these operators are of lower order than what a simple counting of orders (with 
-1 for the order of Vo) would give: 

Lemma 4.7 Let T be smooth. Then the operators curlr VoVr and divr Vo curlr are 

linear and continuous from iT*(r) into itself. 

PROOF. These results are due to the equalities (|3.7|) . One can write (see |28| page 240]): 
curlr VoVtu(x) = J n(x) ■ curP {(7(0, \x — y\)Vru(y)} da(y) 

{(n(x) - n(y)) x V*G(0, \x - y\)} ■ V T u{y)da{y) 
— Vo curlr Vr«- 

The second term on the right hand side vanishes, and the kernel 

(n(x) - n(y)) x V x G(0, \x - y\) 

has the same weak singularity as the fundamental solution G(0, \x — y\). We deduce the 
lemma using similar arguments for the other operator. ■ 

As a consequence, the operators C\\ and C22 are of order -2, the operators C12 and 
£21,1 are of order -1 and the operator 621,2 is of order 1. Therefore, Cq is a compact 
perturbation of 

C12 
C 2 i,a 

By definition of C Kj o, the operator C Ki o can be written as: 



22 



( "'" 1 -«C21,l + « ^21,2 ~kC 2 2 

3 1 

The second term on the right hand side is compact on H2(T)/M. x i?2(r)/R. 
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Since C K) g is a compact perturbation of 
the sum C Ki $ + pC Ket o is a compact perturbation of 

/ -(/Cj + /JKe) 
V (KT l + p^- 1 ) 

The operator C Ki fiC Kei Q is a compact perturbation of 

/ -KiK' 1 \ 2 

Remark 4.8 Notice that lemma \4^\ is n °t true in the Lipschitz case. Nevertheless, one 
can use the Helmholtz decomposition also for a Lipschitz boundary. One only has to 
replace the space iJa(r) by the more general space 

H(T) = {ue if 1 (r), A r u € H~^(T)} , 

see [H [7] . Then a large part of the previous arguments is still valid. For example, the 
operators C\\ and C22, being of order — 1, are still compact from 7-L{T) and ff2(T), re- 
spectively, to themselves. By the compactness of the embedding rL(T) H2 (T) we deduce 
that the operator (721,1 is still compact from rl(T) to Hz(r"). The complete proof of the 
theorem does not go through, however, because it uses the compactness of M K . 

Lemma 4.9 For smooth V , the operator Cg 2 is a compact perturbation of — ll. 

PROOF. It suffices to consider the principal part of (|3.9p , ■ 

Collecting all the results, we find that S\ is a compact perturbation of 

/ \b(p + KiK~ l ) - \a(Ki + pK e ) 

V \b{p + kT 1 K e ) + \a(nr l + pK- 1 ) 

We recall that p = ^ eKl . The matrix written above is invertible if: 

PiK e 




h{p + KiK~ x ) - ^a(Ki + P K e ) + o l(6-2aK e ) fl + ^ ^0 

and 

lb(p + K 7 1 K e ) + ~a(K7 1 + p K J 1 )^0 1(6^ + 20)^1 + ^^0. 

Since the operator Cg is invertible, we conclude that under the conditions of the theorem 
the operator Si is Fredholm of index zero and therefore S too. The theorem is proved. ■ 
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Using similar arguments we obtain the following theorem. 
Theorem 4.10 Assume that 

(i) the boundary T is smooth and simply connected, 
(ii) the constants a, b, p e , pi, K e and m satisfy: 

(bK e ~ 2a) ^ 0, (l + — ^ ^ 0, (b + 2aK e ) ^ and (l + ^-\ \ ^ 0. 

V H) V IHKiJ 

_i 

Then T is a Fredholm operator of index zero in H x 2 (divr,r). 
PROOF. The operator T is a compact perturbation of 

Ti = a + pC Kii0 C Ke A - h - (C Kei0 + pC Kt , ) C *. 

_ i 

We use again the Helmholtz decomposition of H x 2 (divr, T) and rewrite the operators C Ki o 

3 1 

as operators defined on H2(T)/M. x i?5(r)/lL Collecting the results from the previews 
proof we found that the term (C Ke ,o + pC-Ki,o)C$ is a compact perturbation of 

If (n e + pm) 

Finally the operator Ti is a compact perturbation of 

if a (l + pKiK^ 1 ) - \b (n e + pKi) 
4 V a (l + pK e n7 X ) + \b (k' 1 + pnr 1 ) 

Note that under standard hypotheses on the materials and for real frequencies, the ma- 

2 

terial factors such as (l + — ) and (l + ^ eK ^ ~) are always non-zero. 

Pi Pi^e 

Remark 4.11 Thanks to the explicit representations (|4.23p and (|4.24p one can deduce 

_ i 

Garding inequalities (positivity modulo compact perturbations) in H x 2 (divr,r) for S (via 
the bilinear form B) and for T (via the {^-duality pairing) in the case of a domain dif- 
feomorphic to a ball. 

When r is a only a Lipschitz continuous surface, one can still prove that the operators 
S and T are Fredholm operators of index zero, if one imposes more restrictive hypotheses 
on the physical parameters. We have the following Garding inequalities. 

Theorem 4.12 Assume that 
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(i) fi e , Hi, K e and K{ are positive real numbers. 

(ii) a = 1 and b = —in with n E R, n > 0, 

Then the operators S and T satisfy the following Garding inequalities 

Im(£(SjJ) + C5 (j,j)) >C 5 ||j|| i 

H x 3 (div r ,r) 

Re(£(T<7*-HJ) + CTaj)) >Qr||j|| i 

H x 2 (div r ,r) 

where cs and ct are compact bilinear forms and Cs and Ct are positive real constants. 

PROOF. According to the definitions (|4,3p - (|4.7p and Lemma |3.7[ the operator S is a 
compact perturbation of 

Si = -irj (p( - |I + M ) 2 +C Ki , C Ke ,o) C * - § ( P C Kefi + C Kifi ) + P C Kefi M + M C Ki>0 . 
_i 

Let j E H x 2 (divp, T). The term B({pC Ke fi + C Kii o)j,j) is real. We also have 

B(C Ke , M j,j) = -B{M ],C K Ji) = B(i,M C Kefi ]) = -B(M C Kefi H) 
The term MqC Kc q is a compact perturbation of M Ke C Ke and M Ke C Ke = —C Kf ,M Ke , thus 



B(C Kefi M j,j) = ^(C Kej0 M j,j) + compact = Re^(C ttei0 M j,j) + compact 

In the same way we have 

B(M Q C Kit0 J,j) = Re B(M C Kifi l j) + compact. 

Using now the Calderon relations (|3.9p and the fact that Cq = iC K fi for k = i, we see 
that MoCq is a compact perturbation of — CqMq. It follows that ( — |l + Mo) 2 Cq is a 
compact perturbation of ( — |l + Mo)Cq ( — ^ — Mo) and we have 

B (( - ±1 + M )C *( " \l - M ) jj) = B (C *( - ±1 - M ) j, ( - ±1 - M )j) 
Since Cq is elliptic for the bilinear form B we have with Lemma [3.81 

-B (( - |I + M ) 2 C * j, j) + dOJ) = —B (C *(±I + M ) j, (±1 + M )j) 

> ai||(il + M )j|| 2 1 

H x ^(div r ,r) 

> «2 1 1 j || 2 _ 1 -C 2 (j,j), 

H x ^(div r ,r) 

where ci(-, •) and C2(-, •) are compact bilinear forms and ct\ and 02 are positive constants. 

Now for brevity write Sq = n x Vo and Tq* = curlr Vo divr . Taking into account 

_ 1 _i 

that (T *) 2 = and that (S%) 2 : H x 2 (div r ,r) H x 2 (div r ,r) is compact (it maps 
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continuously into H x (divr,r) which is compactly imbedded in H x 2 (divr, T)), we get 
from the definitions 

C Ki fiC Ke)0 Co = -ihk^S^TqSo - k^KsTqSqTq + compact. 

The terms SqTqSq and Tq SqTq give positive contributions, namely we have 

-B{S* T£S* j, j) = ^(curlr V ]) ■ V (cnvl r V j) > 

-B(r *5 *T *j,J) = ^(curlr V div r J) • Vb(curl r V div r j) > 
Therefore there exists a compact bilinear form C3 such that 

-BfC^oC^oCS j, j) + c 3 (j, j) > . 
Collecting all the results, we can write 

Im/3(Sj,J) = V (-B (C *(|l + M ) j, (±1 + Mo)]) - B(C Kii C Kei0 C^ j,J)) - c 5 (j, j) 

where_ -c 5 (j, J) = r/(ci(j,J) + c 2 (j,J) + c 3 (jj)) + ImS((pC Kei oM + M C Kl , ) j, j) + S((S - 
Si) j, j) is a compact bilinear form and 

Im(B(Sj,J)+C50M))>^||j|| 1 

H x 2 (div r ,r) 
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with = 7702- Similar arguments can be used for the operator TC, 

5 Integral equations 2 

The second method is based on a layer ansatz for the interior field: We assume that 

the interior electric field E* can be represented either by ^e k 1 or by ^m k . j where the 
_ i 

density j € H x 2 (divp, T) is the unknown function we have to determine. We begin with 
the Stratton-Chu representation formula fj3.5|) in J7 C : 

E'(s) = * es1 c Nkb (E s + E mc ) (x) + *M Ke l c D (E s + E- c ) (x) xen c (5.1) 

We then apply the exterior traces 7^ and 7^ and use on both sides of (|5.ip the 
transmission conditions. The result is a relation between the traces of E ! on T: 

lD e-Y D V nc = -pC^N^e + {-\l + M Ke ) lD E\ (5.2) 
p 7 ^E*- 7 ^ e E mc = -C KelD V + p{-\l + M Ke ) lNK E\ (5.3) 

In the scalar case, to construct the integral equations one would simply take a linear 
combination of (|5,2p and (|5,3p . Here we multiply (|5.2p by a and (|5.3p by fcCg and subtract 
to obtain: 

/9L' e 7JV K . E' - N' elD E i = h sur T (5.4) 

_i 

where the operators L' e and iVg are defined for all j € H x 2 (divr, T) by : 
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L' e \ = {aC Ke -bC* {\l + M Ke )}\, 
N' e \ = {-a{\l + M Ke )+bC^C Ke }l 

and 

h = a 7D E mc - bC* o7NK E inc E H7(div r ,r). (5.5) 

_ i 

If E l is represented by the potential *&e k . applied to a density j E H x 2 (divr,r) 

P(x) = -(* Bis j)(x), xeQ, (5.6) 

we obtain 

7^ = C K . j and 7JV(S E 4 = (±1 + M re .) j on T (5.7) 
Substituting (15 .7p in (|5.4jl . we obtain a first integral equation: 

S'j = {/< (±1 + M Ki ) - N' e C Ki ] j = /» on r (5.8) 

_ 1 

This is an integral equation for the unknown j E H x 2 (divr,r). Having solved this 
equation, we construct E l and E s by the representations (|5.6[) in Q and 

ES = P (*E Ke (|I + M Kj ) j) (x) + (^M Ke C Ki j) (x) x E n c . (5.9) 

_ 1 

If E* is represented by the potential ^m k . applied to a density m E H x 2 (divr,r) 

E\x) = -(*M Ki m)(x), xEft, (5.10) 

we obtain: 

7D E l = (i/ + M Ki ) m and 7A r K . P = C Ri m on T. (5.11) 

Substituting (|5.1ip in (|5.4p . we obtain a second integral equation: 

T'm = { P L' e C Kt - N' e (±1 + M Ki ) } m = h on T. (5.12) 

_i 

This is an integral equation for the unknown m E H x 2 (divr, T). Having solved this 
equation, we construct E* and E s by the representations (|5,10p in f2 and: 

E'(x) = p ^>E K C Ki m) (x) + (* MKe (±1 + M Ki ) m ) (x), x E n c . (5.13) 

Contrary to the preceding method from the previous section, the two integral equa- 
tions are not equivalent, in general. The following theorem corresponds to Theorem 14.11 
The proof is similar to the scalar case. 

Theorem 5.1 We assume that K^ e is not an eigenvalue of the associated interior problem 

(BHD. 

_ 1 

If] E H x 2 (divr,r) solves (|5.8p . then E* and E s , given by (|5.6p and (|5.9p respectively, 

solve the dielectric scattering problem. 

_i 

Ifm E H x 2 (divr,r) solves (|5.12p . then E l and E s , given by (|5.10p and (|5.13p respectively, 
solve the dielectric scattering problem. 
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PROOF. In each case the integral representations of E l and E s satisfy the Maxwell 
equations and the Silver-Miiller condition. It remains to prove that the transmission 
conditions are satisfied. We prove it for the equation (|5.12p . the arguments being similar 
for (PES) . 

Assume that m solves (I5.12p which we rewrite as : 

a { P C Ke C Kt m + (±1 + M Ke ) (±1 + M Ki ) m - lD E mc ) 

(5.14) 

-bC* {p [\1 + M Ke ) C Kt m + C Ke (II + M Ki ) m - 7 7V Ke E mc } = 0. 

Then, using the integral representation (|5.13p of E s , we obtain : 

( 7 |,E S + 7 f,E mc - 1D V) = - P C Ke C Ki m - (±1 + M Ke ) (±1 + M Ki ) m + 7D E^ C , 
(7k e E s +7k e Em - ^ E*) = -p (il+M Ke ) C Ki m - C Ke (|l+M Ki ) m + 7 ^E mc . 

We have to show that the right hand sides of these equalities vanish. 
We introduce the function v defined on by 



v(x) = - P ^! EK C Ki m - v?M Ke (|I + M Ki ) 



m 



By equation (|5.14p we have a 7 £>v — frC^Ar^v = 0. Since E mc satisfies the Maxwell 
system curl curl v — k^v = in Q, v satisfies it, too. By hypothesis, Kg is not an 
eigenvalue of the associated interior problem, which implies v = in O. In particular, 
7£>v and 77v Ke v vanish, which shows that the above right hand sides are indeed zero and 
that the transmission conditions are satisfied. ■ 

Theorem 5.2 Assume that the hypotheses of Theorem \2.1\ are satisfied and that Kg is 
not an eigenvalue of the associated interior problem (|4.17p . Then the operators S' and T' 
are infective. 

PROOF. We prove the result for the operator T', similar arguments being valid for S'. 

_i 

Assume that mo € H x 2 (divr,r) solves the homogeneous equation : 

T'mo = pL' e C Ki m Q - N' e (±1 + M Ki ) m = 0. (5.15) 

We want to show that mo = 0. 
We construct Vi and v 2 as follows: 

v 2 (x) = p(y EKe C Ki m )(x) + ($M Ke {±I + M K Jm ) (x), x £ fl c , 

and 

vi(x) = -(^M K .m Q )(x), xeQ. 

By Theorem 15. 1| these functions solve the homogeneous scattering problem (i.e. when 
pnc = anc [ therefore Vi = in and v 2 = in £l c . Now we define 

v(x) = -($M K m )(x) x£fL c 
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We have 7^ v = C Ki mo = 7jv k . v 1 = 0. Since v satisfies the Silver-Miiller condition, we 
have v = in U c . Thus v = is R 3 and [7d]v = no = 0. ■ 

Remark 5.3 The operators S' and T' are the dual operators of S and T, respectively, for 
the bilinear form B. Therefore they are Fredholm of index zero under the same hypotheses 
as those given in Theorems \4-5\ and \4-12 



In order that each of the four integral equations admit a unique solution for all positive 
real values of K e , we will now give an example of how to choose the constants a and b 
such that the associated interior problem does not admit any real eigenvalue. 

We summarize all the previous results in the final theorem. 

Theorem 5.4 Assume that the boundary T is smooth and simply connected and 
(i) K e is a positive real number or Im(K e ) > and Re(Kg) 7^ 0, 



(ii) a = 1 and b 



ir) with r) <E R\{0} i/^el 

—ir\K e • sign(Im(Kg)) with rj € R, rj > otherwise, 



2 

(ill) — -I, 1 ^ -1. 

_l 

Then the operators S, T, S and T' are invertible in H x 2 (divr,r). Moreover, given 
the electric incident field E mc € Hi oc (curl, M 3 ) ; i/ie /our integral equations (|4.7I) . (|4.9p . 
(|5.8p . (|5.12p eac/i /jawe a unique solution, and the integral representations { (|4.ip . (|4.1ip |. 
{ (I5.6p . (|5.9p | and { (|5.10p . (|5.13p | o/E* and E s giue the solution of the dielectric scattering 
problem. 

If r is onZy Lipschitz, then the conclusions remain valid if the conditions (i) to (Hi) are 
replaced by the more restrictive assumptions 

(iv) fj, e , Hi, K e and Ki are positive real numbers. 

(v) a = 1 and b = —if) with rj £ M, r] > 0, 



6 Discussion 

In this paper we have described and analyzed modified boundary integral equations to 
solve a radiation problem for the Maxwell system that are stable for all wave numbers. 
Generalizing the approach of Kleinman and Martin to the Maxwell system by employing 
a suitable regularizing operator introduced by Steinbach and Windisch, in Section H] we 
have derived two boundary integral equations using an ansatz for the exterior field and 
in Section [5] we have derived two integral equations using an ansatz for the interior field. 
Note that if it is only the exterior field that is of interest, one can choose an integral 
equation which gives a simple representation for E s , e.g. (|4.7|) or (|4.9p . This choice was 
used in the PhD thesis |21| for an application to a shape optimization problem involving 
the far field pattern |12j . For numerical results using this method, we refer to [21 J. 
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In [22j P. A. Martin and P. Ola established the existence and the uniqueness of the 
solution to an integral equation analogous to (|4.7p for all real values of the exterior 
wave number by adapting a regularization method that was introduced by Kress [9] 
in the framework of spaces of continuous functions, namely by using the operator j h-> 

n x Vq j m the place of our Cg . This technique would not yield four families of Fredholm 

_ i 

boundary integral operators of index zero in H x 2 (divr,r) since the invertibility of the 
regularizing operator is needed in our arguments (see the proof of Theorem 14, 5p . A 
more interesting advantage of the operator Cg is the possibility to use our regularization 
method on Lipschitz boundaries since this operator still is elliptic. 

Numerical analysis using similar CFIE-based methods for the scattering of homo- 
geneous penetrable objects are presented in |33| . The proposed integral equation also 
contains double and triple operator products. Stable discretization of these operator 
products can be obtained by multiplying matrices arising from the discretization of the 
various operators using specific basis and testing functions [lj. As is shown there, pre- 
conditioners from the same class of operators are also easily constructed. The conclusion 
is that appropriately preconditioned CFIE-based single-source formulations are more ef- 
ficient than the coupled integral equations. 
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